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Limits on parameters describing physics beyond the Standard Model are presented. The most 
general Lorentz invariant effective Hamiltonian at the quark-lepton level involving vector, scalar 
and tensor operators has been used. The fits have been done using the most precise and up to date 
experimental data for correlation coefficients measured in free neutron beta decay as well as the 
Fierz term measured in superallowed Fermi decays. 
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I. INTRODUCTION 

The neutron decay, as the simplest beta decay, well 
described in the framework of the Standard Model (SM), 
is used as a tool to find limits on physics beyond the SM. 
There were many approaches to find limits on parameters 
describing physics beyond the SM, based on observables 
measured in the neutron and the nuclear beta decays, 
with the most recent and general ones included e.g. in 
Refs. [1, 2]. We decided to re-run the analysis for several 
reasons. At first, most of the previous studies were based 
on an effective Hamiltonian at the nucleon-lepton level, 
which makes the distinction between New Physics and 
nucleon structure not so clear. Secondly, we aim to find 
bounds on the parameters which in future will be partic- 
ularly useful in the description of neutrino states in beta 
beam oscillation experiments. Moreover, our new fits use 
the most precise and up to date experimental data from 
correlation coefficients a, A, B measured in free neutron 
beta decay (see TABLE I) and the Fierz term measured 
in superallowed Fermi decays (the bp, see Ref. [3]). This 
work continues and extends our previous researches from 
Ref. [4]. 

Similarly as in Ref. [5] we would like to present an 
analysis starting at the quark-lepton level, with a par- 
ticularly simple distinction between the New Physics pa- 
rameters and the nucleon form factors in the vector part 
of the interaction. However, in contrary to what was done 
in Ref. [ ], we performed a self-consistent least squares 
analysis to the free neutron decay data, putting careful 
attention to the still unknown value of the Fierz term 
b. Our paper is divided, except for this short introduc- 
tion, into five sections. In Section II we give a short 
introduction into the beta decay formalism needed for 
our purposes. In Section III we show the definitions of 
the fitted parameters and in Section IV the details of the 
data analysis are given. In Section V we present our re- 
sults and the conclusions are summarized in Section VI. 
In the Appendix, the formulas for decay parameters in 
the free neutron decay are presented. 



II. FORMALISM OF THE NEUTRON BETA 
DECAY 

The most general, derivative-free and Lorentz invari- 
ant four fermion contact interaction for neutron decay 
at nucleon-lepton level was first introduced by Lee and 
Young [6] and later used in the canonical work by Jack- 
son, Treiman and Wyld [7] to obtain formulas for com- 
plete set of decay parameters. Such general interac- 
tion involving scalar, vector and tensor terms can be 
parametrized in many ways (not only as in Refs. [6, 7]). 
Herein, we will use the parametrization introduced in 
Ref. [5] with the interaction Hamiltonian at the quark- 
lepton level in form 

n , 

n e = 4 12 12 {MiU^PmuPid 

i=l k,l=LM. ^ 

+a k iU^e-f tl P k iy l u^Pid 
+a k kU^e^P k u i u^P k d\+R.c, (1) 

where P L = \ (1 - 75), P R = \ (1 + 75), = § [7^ 7„] 
(the metric and gamma matrices are the same as e.g. in 
Ref. [8]), u and d are quark fields, e stands for the electron 
field and Vi is the i-th neutrino field with a certain mass, 
U L and U R are 3 x n mixing matrices 1 for the Left and 
Right-handed neutrinos respectively (we work in the ba- 
sis in which mass matrix of charged leptons is diagonal) . 
We assume in our work that a k i, A k i, a k k for k,l — L, R 
are real, contrary to Ref. [•">]. The SM is restored when 
n = 3 (then U L is the well known Pontecorvo-Maki- 
Nakagawa-Sakata matrix), a k i = 0, A k i = 0, a k k = 
for k,l — L,R except a^L = o-ll — KtdGWv2, where 
Gf is the usual Fermi constant and V u d is the element of 
quark Cabibbo-Kobayashi-Maskawa mixing matrix. 
When calculating amplitudes for the neutron beta de- 



1 Comparing the notation used in Ref. [5] we substitute U with 
U L and V with U R . 
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cay we used the following relations [ ] 

gvUpj^Un = (p\ujpd\n) , 
gAUpj^Un = (p|M7 M 7 5 d|n), 
g s u p u n = (p\ud\n), 
gTUp(Jp V u n = (plua^dln), 



(2a) 
(2b) 
(2c) 
(2d) 



where gi = gi (q 2 w 0) for i — V, A, S, T are the form fac- 
tors at small four-momentum transfer q = p p —p n and (p|, 
p p , u p (|n), p n , u n ) are the proton (neutron) state, four- 
momentum and the ordinary free-particle Dirac wave 
function. From conserved vector current hypothesis, ne- 
glecting small corrections, we know that g v — 1 (see 
Ref. [ ] for a review). There were many attempts to 
calculate gA and gr in lattice QCD, which were summa- 
rized and supplemented with the first calculation of gg 
in Ref. [9] as follows: (i) the central value of gA is in 
the range 1.12 < gA < 1-26 and when we include the 
errors (that are not competitive with the experimental 
uncertainty on the gA/gv ratio — see section V) we can 
conclude that 17,4 can be roughly in the range from 1.1 
to 1.34, (ii) the values of form factors beyond SM are 
g s = 0.8 ± 0.4 and g T = 1.05 ± 0.35. In the later deriva- 
tions in this work gs, gr, gv an d gA are treated as free 
real parameters (we assume time reversal invariance of 
strong interactions). 



over kinematically allowed anti-neutrino states) can be 
accessed through measurements of the total decay width 
- it is however beyond our present interest. The D 
correlation coefficient is mentioned here only for com- 
pleteness, since D = upon approximations under con- 
siderations and for real parameters (experimentally one 
has [10]: D = (-1.2 ± 2.0) x 10~ 4 ). The B correlation 
coefficient has the form of 



B = B + b 



v E P ' 



(4) 



where the formulas for Bq and b v as well as for the cor- 
relation coefficients a, b, A and the common factor £ are 
listed in the Appendix — they are expressed in terms 
of the following parameters (similar as in Ref. [ ] ) for 
k,l = L,R 



A = 
Vrl = 
V RR = 

Sm = 



gA &LL — CLLR 
gv a L L + CLLR 

a-LlARL — O'LRflRR. 
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in 
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gs 
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kl 



gv O-LL + CLLR 



(5a) 
(5b) 
(5c) 

(5d) 



III. DECAY PARAMETERS 

Using Eq. (1) we calculated the five- fold differential de- 
cay width for polarized neutron in its rest frame without 
measurement of final electron and proton polarizations. 
As in Ref. [7] , in our calculations we have neglected QED 
corrections (unless otherwise stated) and recoil effects 
(for a review see e.g. Ref. [9]), obtaining 



dT 



dE„dVt„dVL. 
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(3) 



where Sl e and f2„ denote the solid angles of electron and 
anti-neutrino emission, m e , p e , E e are, respectively, the 
mass, momentum (p e = \p e \) and total energy of an elec- 
tron, p v and |p„| = E v — E — E e are the anti- neutrino 



momentum and energy, Eq 



i p (that is a dif- 



ference between neutron and proton masses) is the maxi- 
mum value of E e and A„ is the neutron polarization vec- 
tor. 

The value of Gp = 2(a LL + a LR f g^'i P^ 2 (the 
prime in the sum indicates that the summation runs only 



where 
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9t 
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gv a L L + a LR 



A~ki = A 



k 12 
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Otkk = Ctkk\/J2i \ U t 



Tk 12 



(5e) 

(6a) 
(6b) 
(6c) 



with summation running only over kinematically allowed 
anti-neutrino states. For a^R = it is easy to sepa- 
rate the dependency of decay parameters on the ratio of 
nucleon form factors A = gA/gv from the rest of param- 
eters. 

For scalar couplings, in general, the formulas for cor- 
relation coefficients (see the Appendix) depend only on 
two combinations 



sl — Sll + Slr, 

Sr = Srr + Srl- 



(7a) 
(7b) 



Note that setting clrr = results in V R k = aRkja^L 
for k = L,R and our previous parametrization [ ] is re- 
stored, except for the old Vlr = <ilr/o.ll parameter 
which no longer appears as the cllr explicitly enters our 



2 The effect of nonzero neutrino masses enters only through pres- 
ence of mixing matrices U L and U R . 



3 Note that Vrr + Vrl = {a RR + a R h)/{aLL + o, L r) and Vrr - 
Vrl = (3.RR - a RL )/(a LL - a LR ). 
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present parametrization — see Eqs. (5) (note also the 
small change of notation U — U L and V — U R ). In fact, 
even for a LR ^ the results presented in plots in Ref. [4] 
can be interpreted in terms of our present parametriza- 
tion provided that gA/gv m those plots is replaced with 
A, as defined in Eq. (5a). 



IV. LEAST SQUARES ANALYSIS 

Following the approach applied in Refs. [1, 2, 11], the 
measurements of a, A, B are interpreted as measure- 
ments of the respective quantities 



PAR. VALUE ERROR (W" 1 } PAPER ID 



adW- 1 )) 

Miw- 1 )) 

mw- 1 )) 



1 



A 



l + biW- 1 }' 
1 + biW- 1 ) 



(8a) 
(8b) 
(8c) 



where (W _1 ) — m e (i? ( r 1 ). This procedure arises mainly 
from the fact that experimentalists analyze their data 
assuming 6=0 and 6„ = 0. Note that in the SM and for 
some combinations of parameters of physics beyond SM 
this assumption is valid. 

The % 2 , which is minimized with the fit procedure, is 
of the form 



E 

i 

E 

3 

E 



a((W- 



A, 



Sa,i 

AHW- 1 ) 



SA j 
5B k 



(9) 



where the selected data are presented in TABLE I: a^, 
Aj, Bk and 8di, SAj, 6B k denote the central value and 
the error of the respective decay parameter in a certain 
experiment. We follow the PDG [10] data selection, but 
(i) we have used the corrected value for Ref. [19] given in 
Ref. [17], (ii) we added new measurements of A [17, 18] 
and dropped older measurements of this decay parameter 
j 3] since they are poorly consistent with the newer 
ones and (iii) we kept only the most precise measure- 
ments of a and B (Joi/a, < 6% and SBk/B^ < 2%). 
When both statistical and systematic errors were re- 
ported separately, these two errors were added in quadra- 
ture. In the case of asymmetric errors the larger of the 
reported errors was taken. 

Part of the values of (II 7 " 1 )^ = m e (E~ 1 ) i , present in 
TABLE I, were taken from Ref. [1]. The remaining ones 
were calculated using 



dE, 
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(10) 
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-0.11996 


0.00058 
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MUND 


12 [18] 




-0.11942 


0.00166 


0.557 


LIU 


10 [17, 19] 




-0.1189 


0.0007 


0.534 


ABELE 


02 [20] 


B 


0.980 


0.005 


0.599 


SCHUMANN 


07 [24] 




0.967 


0.012 


0.600 


KREUZ 


05 [25] 




0.9801 


0.0046 


0.594 


SEREBROV 


98 [26] 




0.9894 


0.0083 


0.554 


KUZNETSOV 


95 27 



TABLE I. The experimental values of correlation coefficients 
measured in free neutron beta decay. All PAPER ID names 
were taken from PDG [10] except the new ones for Ref. [17] 
and Ref. [18]. 



where E™" %n and E™ ax are the limits of the energy 
accepted in the experiment. The differential decay 
width was calculated within SM, therefore (compare e.g. 
Ref. [14]) 



d^SM 

dE, 



= (9v + 39a) 



G%\V ud \< 



2tt 3 



Pe E e (E -E e ) 2 F(E e ), 
(11) 

where the Fermi function F(E e ), that is the leading order 
QED correction, was approximated (for an exact calcu- 
lation see Ref. [12]) as in e.g. Ref. [13] to the form of 



F(E e 



2naE e /p e 



1 



-2-KaE e /p e 



(12) 



The Vrl, Vrr, sr, Trr enter quadratically or as mixed 
terms between pairs of these couplings in the formulas for 
the correlation coefficients (see the Appendix), whereas 
Sl and Tll enter also linearly. Therefore, the x 2 function 
in Eq. (9) reveals the symmetry 

X 2 (\s L ,T LL ,V RL ,V RR ,s R ,T RR ) = 
X 2 (X, s L ,T LL , -Vrl, -Vrr, Sr, -Trr), (13) 

The linear terms appear only in b and b v so these quan- 
tities are of special importance from the point of view of 
setting limits on sl and Tll- 



V. RESULTS 

a. The A value. Let us first consider the case when 
all aki, Am, ctkk parameters are zero except cill and <1lr- 
This results in Vrl = 0, Vrr = as well as T kk = and 
Sfc = for k — L, R so that the only nonzero parameter 
is A given in Eq. (5a). In this case b = as well as b v = 
and the formulas for decay parameters simplify to the 
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well known SM expressions 



A = 



B = 



A — 
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(14a) 
(14b) 
(14c) 



whereas A can go beyond its SM form of A = qa/ gv (see 
Eq. (5a), note that in our convention A > contrary to 
e.g. PDG [ ]). In this case, the one-parameter fit to 
the data presented in the TABLE I is performed, which 
results in x 2 min — 9.542 (the value of x 2 at minimum) 
with 



A = 1.2755 ± 



0.0011 (68.27% C.L.), 
0.0018 (90% C.L.), 
0.0022 (95.45% C.L.). 



(15) 



The PDG average is [10]: A = 1.2701 ± 0.0025 (the error 
was scaled by PDG by 1.9, we changed the sign to meet 
our convention) which differs from our value because of 
different data selection, mainly of the A decay parameter. 

b. Many-parameter fits. Next we consider cases 
when one or two of Vrl, Vrr, Tkk, Sk parameters are 
nonzero together with A, that goes beyond its SM form 
as before. Results of such two and three-parameter fits 
are presented in FIGs. 1-4. Let us start with a few gen- 
eral remarks: (i) we list only those arguments of \ 2 that 
are actually fitted — the rest are set to while A, if not 
fitted, is set to its central value given in Eq. (15), (ii) the 
two-dimensional plots in the case of three-parameter fits 
are obtained by intersecting the corresponding three- 
dimensional x 2 volume with a plane that includes the 
three-parameter x 2 minimum point and is parallel to the 
respective planes spanned on the main axes in the pa- 
rameter space. Enclosing general remarks we describe 
particular fits in more details below. 

First, we made fits to the neutron data presented in 
the TABLE I in all possible two and three parameter 
combinations in which 6 = and b v = 0. These results 
can be divided into two main groups: (i) the fits when 
the only nonzero parameters are Vrl ° r Vrr or both of 
them simultaneously (vector couplings fits) — see FIG. 1, 
(ii) the fits when the only nonzero parameters are sr or 
Trr or both of them simultaneously (right tensor and 
right scalar couplings fits) — see FIG. 2. Combining 
vector couplings (first group) with right tensor or right 
scalar couplings (second group) results in b and b v being 
not identically zero. In the case of two-parameter fits 
there are 2 equivalent minima because of the x 2 symme- 
try given in Eq. (13). In the case of three-parameter fits 
there are 2x2 = 4 equivalent minima: the minimization 
procedure found two equivalent minima corresponding to 
the different values of A and, for each of these values, we 
have two sets of Vrl, Vrr or sr, Trr parameters from 
X 2 symmetry in Eq. (13). 

Finally, we made fits when the only nonzero parameter 
is sl or Tll as well as both of them simultaneously (left 



tensor and left scalar couplings) using only neutron data 
as before — these results are presented in FIG. 3. There 
is only one minimum in each fit and in all cases b as well 
as b v are not identically zero. Moreover, we report that 
the calculated values of b and b v are rather big (see e.g. 
[28]) especially when all parameters are fitted together 
as it can be seen in the TABLE II — this causes strong 
dependency of the results on the (W~ 1 )i values. To re- 
duce this effect we added to the data set a measurement 
of b in superallowed Fermi decays [3] 



b F = -0.0022 ± 0.0026. 



(16) 



We therefore extended our previous x 2 m Eq. (9) with 
an extra term 



2 (bp + 0.0022)^ 
X + (0.0026) 2 



where the general formula for bp in given by 4 
-2 [s R (V RR + V RL ) 



b F = 



(Vrr + VrlY + sI + sI + 1 



(17) 



(18) 



These results are presented in FIG. 4 and the calculated 
values of b and b v are listed in the TABLE II. From 
Eq. (18) we see that in the case under consideration 
(i.e. when Vrl = 0, Vrr = 0, s R = and T RR = 0) 
bp = — 2sl/(1 + s 2 L ). Thus, including the experimental 
value for bp really affects only fits in which sl is one of 
the nonzero parameters, while x 2 (^^ll) is affected just 
by an overall shift of (0.0022) 2 /(0.0026) 2 « 0.716 with 
respect to the same x 2 calculated without bp. When 
comparing FIG. 4 with plots in FIG. 3, the limits on 
Sl are stronger and the ones on Tll are practically un- 
changed. 

It is difficult to compare our results with the recent 
ones in Refs. [1, 2] because of different parametrizations. 
If we set aLR = we can easier compare our results 
with the analysis in Ref. [ ]. When only one beyond 
SM parameter is fitted together with A = gA/gv to free 
neutron data we conclude that our constrains are not 
as strong as those coming from joint analysis where also 
pion and nuclear beta decays are included. 

We would like to point out that, in cases when A is fit- 
ted together with only one parameter beyond the Stan- 
dard Model, there is no significant difference with respect 
to our previous results [4] except for a shift in A, caused 
by the different selection of A decay parameter measure- 
ments. As expected, the change of the data selection 
influences the result of the one-parameter fit on A. 



4 The formula for bp can easily be derived from formulas obtained 
in Ref. [7], after appropriate change of the parametrization in 
the form given in Ref. [ ] . We assume that in the case of nuclear 
decays, that were studied in Ref. [3], to obtain the quoted bp 
value the same anti-neutrino states are kinematically allowed as 
in the case of the free neutron beta decay. 
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,,2 
X 


without bp 


with bp 


X'{\T LL ) 


-0.022 < 6(0.1003) < 0.19 
— 0.020 < o„(0.092oj <U.1( 


the same as 
without bp 




-0.052 < 6(-0.0276) < 0.0096 


-0.00075 < 6(0.000352) < 0.0015 
n nnnofi <r h (n oooaaq^ <^ n nm q 

u.uuuyu \ Uj/^u.uuu^yj \ u.uuiy 


X 2 (sl,Tll) 


-0.024 < 6(-0.00656) < 0.011 
-0.037 < &„ (-0.0142) < 0.0091 


-0.015 < 6(0.00135) < 0.016 
-0.014 < 6^(0.00137) < 0.015 


X 2 (X,s l ,T ll ) 


-0.28 < 6(0.0821) < 0.58 
-0.28 < 6„(0.0742) < 0.58 


-0.045 < 6(0.1019) < 0.20 
-0.041 < 6„(0.0942) < 0.18 



TABLE II. The maximal and minimal values of b and b v calculated from parameter combinations at 95.45% C.L.. The values 
in parentheses correspond to \ 2 minima. 



In the case of three parameter fits with x 2 (A, Vlr, Vrr) 
and X 2 (^j S R> Trr) we obtained four equivalent minima. 
Two of them give the A which differs from the SM value 
at the level of 3cr, and the next two where the difference 
is much larger. This discrepancy is also large in three 
parameter fits with x 2 (A, sl,Tll). 

VI. CONCLUSIONS 

We found limits on parameters describing physics be- 
yond the Standard Model in the neutron beta decay in- 
cluding the last experimental data. The fits were done 
using the parametrization in which the influence of New 
Physics at the quark-lepton level can be easily separated 
from the nucleon structure in the vector part of the inter- 
action. Our result for one parameter SM fit of A slightly 
differs from that of PDG because of different data selec- 
tion. In the presence of New Physics one can find values 
of A at the \ 2 minimum that are relatively close to that in 
the SM and other which deviate significantly form the SM 
expectations. We confirm what was already mentioned in 
the literature that free neutron decay alone cannot give 
us the strongest limits for exotic couplings. Therefore we 
are looking forward for new, even more precise experi- 
mental data, especially of the still unknown values of b 
and b v coefficients. 
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Appendix: Formulas for correlation coefficients 

Below we present the formulas for the correlation co- 
efficients a, b, A, B (B = B a + b v m e / E e ) from Eqs. (3,4) 
as functions of the parameters defined in Eqs. (5,6,7). 
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FIG. 1. The results of fits of vector couplings to the data presented in the TABLE I. The list of fitted parameters is indicated 
as the arguments of the \ 2 function in the left upper corner of each plot and the value of the \ 2 a t minimum is shown in the 
right upper corner. The cross marks the position of the x 2 minimum. The solid vertical lines mark areas that correspond to 
the 95.45% C.L. interval of the one-parameter fit on A (Eq. (15)). In the left column, results of two-parameter fits are shown. 
In the center and right columns we present slices of the three-dimensional \ 2 volume in two of four equivalent minima of the 
respective three-parameter fit — the minima corresponding to different values of A are in separate columns. The slices in the 
remaining two minima can easily be obtained through appropriate symmetries according to Eq. (13). 
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FIG. 2. The results of fits of right tensor and right scalar couplings to the data presented in the TABLE I. The arrangement 
of plots and description is analogical to that in FIG. 1. 
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FIG. 3. The results of fits of left tensor and left scalar couplings solely to the data presented in the TABLE I. The first row 
groups slices of the three-dimensional \ 2 volume of the respective three- parameter fit, while in the second we present results of 
two-parameter fits. As in FIG. 1 we list fitted parameters as the arguments of the x 2 function, we quote the value of the \ 2 a t 
the minimum and mark by the cross the position of the x 2 minimum. As before, the solid vertical lines mark the 95.45% C.L. 
interval on A (Eq. (15)). 
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FIG. 4. The results of fits of left tensor and left scalar couplings to the data presented in the TABLE I supplemented with the 
measurement of &f- The arrangement of plots and description is analogical to that in FIG. 3. 
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